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ABSTRACT 

We argue that, at least a fraction of the newly discovered population of ultra- faint 
dwarf spheroidal galaxies in the Local Group constitute the fossil relic of a once ubiq- 
uitous population of dwarf galaxies formed before reionization with circular velocities 
smaller than v c ^ ~ 20 km/s. We present several arguments in support of this model. 
The number of luminous Milky Way satellites inferred from observations is larger than 
the estimated number of dark halos in the Galaxy that have, or had in the past, a 
circular velocity > v^ r , as predicted by the "Via Lactea" simulation. This implies that 
some ultra-faint dwarfs are fossils. However, this argument is weakened by recent re- 
sults from the "Aquarius" simulations showing that the number of Galactic dark matter 
satellites is 2.5 larger than previously believed. Secondly, the existence of a population 
of ultra-faint dwarfs was predicted by cosmological simulations in which star formation 
in the first minihalos is reduced - but not suppressed - by radiative feedback. Here, we 
show the statistical properties of the fossil galaxies in those simulations are consistent 
with observations of the new dwarf population and with the number and radial distri- 
bution of Milky Way satellites as a function of their luminosity. Finally, the observed 
Galactocentric distribution of dwarfs is consistent with a fraction of dSphs being fos- 
sils. To make our case more compelling, future work should determine whether stellar 
chemical abundances of simulated "fossils" can reproduce observations and whether the 
tidal scenarios for the formation of Local Group dwarf spheroidals are equally consistent 
with all available observations. 

Subject headings: cosmology: theory — galaxies: formation — stars: formation 



1. Introduction 

Hierarchical formation scenarios predict that most of the galactic halos that formed before 
reionization (at z = 6 — 10) had masses below 10 8 — 10 9 M . Those that survived to the modern 
epoch, if they were able to form stars, would constitute a sub-population of dwarf satellites orbiting 
larger halos. N-body simulations of cold dark matter (CDM) predict a number of dark matter 
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halos arou nd the Milky Way and M31 that is mu ch greater than the number of observed luminous 
satellites (jKlvpin et al.l Il999l ; iMoore et al.l ll999). This may indicate a problem with the CDM 
paradigm or that feedback processes are very efficient in suppressing star formation in the first 
small mass halos, which remain mostly dark. Recent observational and theoretical advances require 
a re-visitation of the missing galactic satellite problem. In addition, cosmological simulations of 
the formation of the first galaxies have shown most previously known dwarf spheroidals (dSphs) 
to have properties consistent with the surviving first galaxies , and predicted the exi stence of an 
undiscovered, lower surface brightness population of dwarfs (IRicotti &: Gnedinl 120051 ) (hereafter, 
RG05). The recent discovery of a population of ultra- faint dwarfs confirms the aforementioned 
theoretical expectation and allows us to test in great detail cosmological simulations of the first 
galaxies, addressing important theoretical questions on feedback in the early universe and on the 
nature of dark matter. 

The formation of the first dwarf galaxies - before reionization - is regulated by complex feed- 
back effects that act on cosmological scales. These self-regulation mechanisms have dramatic effects 
on the number and luminosity of the first, small mass galaxies, and yet, are unimportant for the for- 
mation of galaxies more massive than 10 s M , that may form before and after reionization. Galaxy 
formation in the high redshift universe is peculiar due to (i) the lack of important coolants, such as 
carbon and oxygen, in a gas of primordial composition and (ii) the small typical masses of the first 
dark halos. The gas in halos with circular velocity smaller than 20 km s _1 (roughly corresponding 
to a mass M ^ 10 8 M at the typical redshift of virialization) is heated to temperatures ^ 10, 000 
K during virialization. At this temperature, a gas of primordial composition is unable to cool 
and initiate star formation unless it can form a sufficient amount of H2 (xh 2 ~ 10 -4 ). Although 
molecular hydrogen is easily destroyed by far ultraviolet (FUV) radiation (negative feedback), its 
formation can be promoted by hydro gen ionizing radiati on emitted by massive stars, through the 
formation of (positive feedback) (jHaiman et al.l ll996). 



It is difficult to determine the net effect of radiative feedback on the global star formation 
history of the universe before the redshift of reionization. The effect of a dominant FUV background 
(at energies between 11.34 eV and 13.6 eV) is to destroy H2, the primary coolant at the start of 
galaxy formation. The FUV radiation emitted by the first few Population III stars is sufficient 
to suppress or delay galaxy formation in halos with circular velocities v c < v^ r ~ 20 km s . 
According to this scenario, most halos with masses < 10 8 — 10 9 M do not form stars and remain 
dark. Therefore, the number of pre-reionization fossils in the Local Group would be expected to be 
very small or zero. However, thi s model does not take into account the effect of ionizin g radiation 
and "positive feedback regions" (jRicotti et al.ll200ll ; lAhn et al.ll2006l ; IWhalen et al.1120081), that m ay 
have a dominant role in regulating galaxy formation before reionization ( Ricotti et al. 2Q02al fbh . 
Simulations including these processes show that star formation in the first small mass halos is 
inefficient, mainly due to winds produced by internal UV sources. This produces galaxies that 
are extremely faint and have very low surface brightnesses. However, according to the results of 
our simulations, a large number of these ultra-faint dwarfs (a few hundred galaxies per co- moving 
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Mpc 3 ) form before reionization at z ~ 7 — 10. Hence, according to this model, the Local Volume 
and the Local Group should contain hundreds of ultra- faint dwarf galaxies. 

The small masses of the first minihalos have two other implications. First, the ionizing radiation 
emitted by massive stars can blow out most of the gas be fore SN-driven winds become important, 
further reducing star formation rates (IRicotti et al.ll2008l ). Second, the increase in temperature of 
the intergalactic medium (IGM) to 10, 000 — 20, 000 K due to H I reionization, prevents the gas 
from condensing into newly virialized halos w ith circular velocities smaller than 10 — 20 km s _1 , 
with a critical value that depends on redshift (|Gnedinll2000l ; lOkamoto et al.l l2008). It follows, that 
dwarf galaxies with v c < 10 — 20 km s _1 stop forming stars after reionization. However, the value 
v"' ~ 20 km s™ 1 that we use to define a fossil is primarily motivated by the fundamental differences 
discussed above in cooling and feedback processes that regulate star formation in the early Universe 
and is not the critical value for suppression of star formation due to reionization. Indeed, IRicotti 
(|2009i ) argues that fossils dwarfs can have a late phase of gas accretion and star formation well after 
reionization, at redshift z < 1 — 2. Thus, a complete suppression of star formation after reionization 
(about 12 Gyr ago) is not a defining property of a fossil dwarf. 

Data on the velocity dispersion of the stars, a, the only observational measure of halo mass, 
shows that a typical dSph has a < 20 km s — 1 . However, the circular velocity of the dark halo 
can be much larger than a if typical radius of the stellar spheroid is much smaller than the dark 
halo radius. For a pre-reionization fossil, according to simulations in RG05, on average we measure 
a/v c ~ 0.5 at formation (see Fig. 2 in RG05). 

RG05 compared the statistical properties of simulated pre-reionization galaxies to observations 
of Local Group dwarfs available in 2004-2005. Based on similarities between observed dSphs and 
simulated galaxies formed before reionization, they argued that many dSphs may be "fossils" of 
the first galaxies. R G05 also predicted the exi stence of a yet undetected population of ultra- 
faint dSph galaxies. (jGnedin &; Kravtsovl 120061 ) (hereafter, GK06) used results from RG05 to 
predict the radial distribution of fossils around the Milky Way and their Galactocentric luminosity 
function. They found a good agreement of simulations with observations for the most luminous 
fossils but, again, a lack of observed ultra-faint fossils, mainly in the outer Milky Way halo. An 
ultra-faint population of dSphs has now been found in the Local Group. As we will show in the 
present work, these new dwarfs have properties in agreement with our simulations of pre-reionization 
fossils. This discovery is certainly one of the most exciting developments in understanding galaxy 
formation in the early universe and has drawn renewed attention to "near field cosmology" as a 
tool to understand galaxy formation. The new galaxies have been discovered by data mining the 
SDSS and surv eys of the halo around M31, r esulting in the di s covery of 12 new u l tra-faint Milky 
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Here, we also argue that solely based on the observed number of Milky- Way satellites, at least a 
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few of them must be a pre-reionization fossil, that means: it must have formed before reionization 
in a halo with circular velocity < 20 km/s. However, some ultra-faint dwarfs may not be pre- 
reionization fossils. As we write, several works have been published that seem to show that the 
observed properties of the ultra-faint dwarf population can also be explained in the context of the 
tidal scenario, that assumes that these galaxies formed after reioniz ation in halos that were much 



more massive and had different properties at the time of formation (jPeharrubia et al.ll2008l ). It is 
intriguing that both the tidal scenario and the pre-reionization fossil scenario are able to produce 
a population of ultra-faint dwarfs that follow very similar statistical trends in terms of size, surface 
brightness, mass to light ration and metallicity- luminosity relation. The jury is still out. 

This paper is laid out as follows. In § [21 we collect published data on the new dwarf population 
and, after correcting for completeness of the surveys, we estimate the total number of Local Group 
satellites (which increases from 32 to about 100). Using the results of published N-body simulations, 
we compare the observed number of luminous satellites to the estimated number of dark satellites 
that have or had in the past a circular velocity > v^, using the results of published N-body 
simulations, concluding that some ultra-faint dwarfs must be pre-reionization fossils. In § [3] we 
show that the properties of the new Milky Way and M31 dwarfs are in remarkable agreement with 
the theoretical data on the "fossils" from RG05 and with their Galactocentric distribution around 
the Milky Way calculated in GK06. In § H] we discuss the implications of the new dwarfs on the 
formation of the first galaxies and the missing galactic satellite problem. 



2. Data and Completeness Corrections 



In Table 1, we summarize the observed properties of the new dwarfs. The new Milky Wa y 



satellites were discovered using SDSS Data Release 4 and 5 (jAdelman-McCarthy et al.ll200d . 120071 ). 
When multiple references are available for a dwarf property, we defer to the measurement with the 
smallest error bars. Excepting Bootes I and II, Canes Venatici I and Leo T, where central surface 
brightness measurements were available, the average surface brightness inside the half light radius, 
ri/2) was used: £y = Ly/(27rr^, 2 ). 



Recent surveys of M31 (jMartin et alJ 120061 ; llbata et al.l 120071 ) have covered approximately a 
quarter of the space around the M31 spiral. The survey have found 6 new M31 satellites. If we 
make a simple correction for the covered are of the survey we find that, including the new dwarfs, 
the estimated number of M31 satellites increases from 9 to 33 ± 10. Two new M31 satellites. 



And XII and And XIV, hav e velocities near or above their host's escape velocity (j Chapman et al 



20071 : iMaiewski et al.1 120071 ) . Both galaxies are classified as dwarf spheroidals and show a lack of 
H I gas, and both are likely on their first approach towards a massive halo. Kinematic data is 
not yet available on these two dwarfs to determine whether their circular velocities are below the 
20 km s _1 threshold, however their currently known properties meet the RG05 criteria for fossils. 



In estimating the completeness correction for the number of the Milky Way dwarfs, one should 
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also account for se lection effects from the limiting surface brightness sensitivity of the Sloan of 
~30 mag arcsec -2 ( Koposov et al.ll2008l ). The sensitivity limit is shown as a solid line in Figure [TJ 
Identification of new satellites depends on the visibility of the horizontal branch in the color- 
magnitude diagram, which, for the typical luminosity of the new faint dwarfs (M y ~ —4) drops 
below SDSS detection limits at Galactocentric distances beyond ~ 200 — 250 kpc (jKoposov et al 



2008). Of the new Milky Way dwarfs, only Leo T is well beyond this distance threshold and nine of 
the eleven new Milky Way satellites are within 200 kpc. We make the most conservative estimate, 
by assuming that we have a complete sample of dwarfs within 200 kpc. Additional selection bias 
for the new dwarfs comes primarily from the limits of the SDSS coverage on the sky. To account 
for this, we apply th e zero-th order correction of multiplying the number of new dwarfs by 5.15 
(jTollerud et al.l l2008). This correction assumes an isotropic distribution of satellites when observed 
from the Galactic center. With this simple assumptions we estimate that the number of Milky 
Way satellites with Galactocentric distance < 200 kpc is about 85 =L 14, including the 29 previously 
know satellites. The error estimate is due to shot noise. 



However, bright satel lites of the Milky Way are distributed very anisotropically (jKroupa et al 



2005 



Zentner et al.l 120051 ) . so the assumption of isotropy may not be a good one. In addition, 



the luminous satellites can be radially biased, so the abundance of the faintest satellites within 
50 kpc may not be easily corrected to larger distances without prior knowledge of this bias. And, 
of course, satellites of different luminosity and surface brightness will have differen t completeness 



limits . These selection biases have been considered in detail in a recent paper by I Toller ud et al 



(|2008l ). This study finds that there may be between 300 to 600 luminous satellites within the 
virial radius of the Milky Way. Their estimate for the number of luminous satellites within a 
Galactocentric distance of about 200 kpc is 120, that is slightly larger than our simple (and more 
conservative) estimate. 



2.1. Number of non-fossil satellites in the Milky Way 

In this section, we use the results of published N-body simulations to estimate the number of 
dark halos in the Milky Way that have, or had, a circular velocity v c > 20 km/s. By definition, 
dwarf galaxies formed in these dark halos are not pre-reionization fossils. If we find that the 
number of observed Milky Way satellites exceeds the estimated number of these massive halos we 
must conclude that at least a fraction of the observed Milky Way satellites are pre-reionization 
fossils. GK06 have estimated that pre-reionization fossils may constitute about 1/3 of Milky Way 
dwarfs, based on detailed comparisons between predicted and observed Galactocentric distributions 
of dwarf satellites. 

It is clear that if we simply count the number of dark halos within the Milky Way virial radius 
with v c ^ 20 km/s, their number is much smaller than the current number of observed luminous 
satellites. However, a significant frac tion of dark halos tha t today have v c < 20 km/s were once 



more massive, due to tidal tripping (jKravtsov et al.1 |2004j) . If the stars in these halos survive 
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tidal stripping for as long as the dar k matter, they may ind eed account for a fraction or all of the 
newly discovered ultra- faint dwarfs. iKravtsov et al.l (|2004l ) favor the idea that tidal stripping of 
the dark matter halo does not affect the stellar properties of the dwarf galaxy. Thus, this model 
is qualitatively similar to our model for pre-reionization fossils, save a rescaling of the mass of the 
dark halos hosting the dwarfs. 

However, it is also possible that tidally stripped halos loose their stars more rapidly than they 
loose their dark matter. Th us, they may quickly tra nsform from luminous to dark halos. Such a 
behavior has been found by iPeharrubia et al.1 (|2008l ). According to this scenario, tidally stripped 
dark halos may not account for the observed ultra-faint population. To summarize, if the number 
of dark halos that have or had in the past v c ^ 20 km/s is smaller than the number of luminous 
Milky Way satellites we may conclude that some dwarfs are fossils. Vice versa, if the number is 
larger, we cannot make any conclusive statement about the origin of Milky Way satellites. 

High resolution N-body simulations of the Milky Way system give the number of dark halos 
in the Milky Way as a fu nction of their circular velocity v c at z = 0. The "Via Lactea" simulation 
bv lDiemand et al.l iftOOi ) finds: 



Ndm{> Vc) = N dm ,20 



20km/s 



(1) 



with Ndm t 20 ~ 27.7 and a « 3. However, a recent work by ISpringel et al.l (|2008l ) (the Aquarius 
simulations) finds a factor 2.5 more satellites at any given v c , i.e., Ndm,20 ~ 69 and a « 3.15. 
Although the Aquarius simulations have higher resolution than the Via Lactea simulation, the 
large disagreement between the two works is due to a systematic difference, possibly related to 
the creation of the initial conditions, and it is not due to the improved resolution. Although the 
Aquarius simulation is likely correct, we will provide predictions for both simulations (we have 
become aware of the Aquarium simulation results, that are not published yet, after this work had 
been mostly completed.) 

To dete rmine the impo r tance of tidal mass loss for satellites aro und the Milky Way we use 
results from IKravtsov et al.l (|2004l ). Figure 5 in IKravtsov et al.l ([20041 ) gives the fraction of halos, 
f(v c ), that presently have circular velocity v c , but some time in th e past had a circular velocity 
< y™ ax = 20 km/s, where as v™ ax = m&x(v c (t)). We approximate the IKravtsov et al.l (|2004l ) results 
for f(v c ) with the power law f(v c ) ~ (v c /20 km s -1 )' 3 , with f3 ~ 3.7. We then calculate the number 
of dark halos N dm {v™ ax > 20 km s _1 ) analytically: 



r20 km/s j iu 

N d m(v? ax > 20 km s- 1 ) = N drn (v c (z = 0) > 20 km s" 1 ) + / dv ~^f( 



N, 



dm,20 



1 + 



O 



(3 — a 



(1 - x 



min ) 



2MN dm , 



20 



where x m i n = v m i n /20 km/s, and v m i n £ < 10 km s 1 roughly equals the mean observed 

velocity dispersion of the stars, < cr # >, of ultra-faint dwarf satellites. The rationale for integrating 
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to 

Vmin is that observed satellites are dark matter dominated and cannot be hosted in dark halos 
that have v c < a*, unless <r* is not a tracer for the dark matter content of the halo (e.g., due to 
tidal heating). 

Using the above equation, we find 73 ± 16 and 182 ± 40 halos with Ndm{v™ ax > 20 km s _1 ) 
within R v i r , for the Via Lactea and Aquarius simulations respectively. B oth these numbers ar e 
smaller than the 300 — 1000 luminous Milky Way satellites estimated by I Toller ud et al.l ((2008). 
Taken at face value, these numbers indicate that a fraction of Milky Way satellites are true pre- 
reionization fossils. However, the number of luminous satellites that exist within the Milky Way is 
highly uncertain beyond a distance from the Galactic center of 200 kpc. 



Based on Figures 11 and 12 in lSpringel et al.l ( 20081 ) and Figure 5 in lTollerud et al.1 (2008]), we 
estimate that roughly half of the Milky Way satellites (within the virial radius R v i r ) are < 200 kpc 
from the Galactic center. Therefore, within 200 kpc we estimate 36 ± 8 and 91 ± 20 dark halos with 
Ndm( v ™ ax > 20 km s _1 ) for the Via Lactea and Aquarius simulations respectively. These numbers 
can be compared to our estimated number of luminous satellites with d < 200 kpc (85±14 satellites) 
and to the estimate by I Toller ud et al.l (|2008l ) (120 satellites). Using the Via Lactea simulation, we 
still find that some dwarfs are true pre-reionization fossils but the argument is weak if we use the 
Aquarius simulation results. In Table 12.11 we summarize the counts for dark matter and luminous 
satellites discussed in this section. 

Although there is considerable uncertainty in our estimates, it can be safely concluded that, 
using the results of the Via Lactea simulation, at least a fraction of Milky Way dwarfs are fossils. 
However, this argument does not hold anymore or is weakened due to recent results (not yet 
published in a refereed journal) from the Aquarius simulations, showing a factor 2.5 increase for 
the number of Milky Way dark matter satellites in any mass range. 



Table 1. Number of observed satellites versus number of dark halos with v™ ax (t) > 20 km/s 
(i.e., non pre-reionization fossils) for the Milky Way 





number of luminous dwarfs 


number of dark halos N^m 


with v™ ax > 20 km/s 


distance from 






Via Lactea sim. 




Aquarius sim. 


Galactic center 


this work 


Tollerud et al. 


today any time 


today 


any time 


d < 200 kpc 


85 ± 14 


120 


14 36 ± 8 


34 


91 ±20 


d < R vir ~ 400 kpc 




from 300 to 600 


28 73 ± 16 


69 


182 ± 40 
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2.2. Peculiar ultra-faint dwarfs 

Almost all newly discovered dwarfs are dSphs with a dominant old population of stars and 
virtually no gas, which makes them viable candidates for being pre-reionization fossils. However, 
there are two notable exceptions that we discuss below. 



2.2.1. Leo T 



With t he gas and young stars of a typical dlrr and the radius, magnitude, mass and metallicity 
of a dSph (llrwin et al.l 120071 : ISimon &: Gehal 1200711. Leo T presents a puzzle. Leo T has a stellar 



velocity dispersion of OheoT = 7.5 ± 1.6km s 1 ( Simon fc Gehal 120071 ) , or an estimated dynamical 



mass of 10 7 M within the stellar spheroid (although its total halo mass may be much larger) 



Leo T shows no sign of recent tidal disruption by either the Milky Way or M31 (|de Jong et al 



2008) and is located in the outskirts of the Milky Way at a Galactocentric distance of 400 kpc. 
Leo T photometric properties are identical to those of pre-reionization fossils. However, if we 
assume that Leo T is a pre-reionization fossil, it is not expected to retain significant gas or form 
stars after reionization. 

How did Le o T ke ep its H I and how did its < 9 Gyr old stellar population form? Work 
by IStinson et al.l (|2007l ) suggests cyclic heating and re-cooling of gas in dwarf halos can produce 
episodic bursts of star formation separated by periods of inactivity. The lowest simulated dwarf 
to form stars has a = 7.4km s _1 , similar to Leo T. However, the increasing of the IGM Jeans 
mass after reionization should prevent gas from condensing back into halos with circular velocity 
v c < 20 km s _1 . Another proposal (Ricotti 2008) is that, although the mass of Leo T at formation 
was ~ 10 7 — 10 8 M (i.e., a fossil), as indicated by its stellar velocity dispersion, its present dark 
matter mass and the halo concentration has increased after virialization by roughly a factor of 
10/(1 + z). This is expected if Leo T has evolved in isolation after virialization, as it seems to 
be indicated by its large distance from the Milky Way. In this scenario, Leo T stopped forming 
stars after reionization, but it was able to start accreting gas again from the IGM very recently (at 
z 1 — 2). This can explain the < 9 Gyr old stellar population and the similarity of Leo T to the 
other pre-reionization fossils. 



2.2.2. Willman I 



Assuming virial equilibrium Wilman I has a dynamical mass within the largest stellar orbit 
(r ~ 100 pc) of 5 x 1O 5 M and a M/L of ~ 470, similar to other ultra-faint dwarfs. However, 
Willman I has a distance fro m the Galactic center, luminosity, surface brightness and size typical 
of a galactic globular clusters (IWillman et al.ll2005al ). The combination of dark matter domination 
with a small core radius, and a mass an order of magnitude smaller than the other SDSS dwarfs 
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suggests Willman I exists in the region between dwarf spheroidals and globular clusters. For our 
discussion of the missing galactic satellites, we treat Willman I as a dwarf because of its large dark 
matter fraction. However, its properties do not agree with being a "fossil" of the first galaxies (see 
Figure Q] in the next section) . 



3. Comparison with Theory 

In this section, we compare the properties of the new dwarf galaxies discovered in the Local 
Group to the theoretical predictions of simulations of primordial galaxies formed before reionization. 
The argument that justifies this comparison is as follows. 

After reionization, due to IGM reheating, the formation of galaxies smaller than 20 km s _1 
is inhibited bec ause the therm al pressure of the IGM becomes larger than the halo gravitational 



potentials (e.g.. iGnedinl 12000 ) . Galaxies formed before reionization stop forming stars due to the 
progressive photo-evaporation of their interstellar medium by the ionizing radiation background 
(although most of the ISM was already lost due to UV driven galactic winds and SN explosions). 
Hence, if pre-reionization dwarfs do not grow above v c = 20 km s _1 by mergers, thei r stellar popula- 



tion evolves passively as calculated by stellar evolution models such as Starburst99 (jLeitherer et al 



19991 ) . We define such galaxies as pre-reionization "fossils". 



Clearly, we do not expect two perfectly distinct populations of fossils with v c < 20 km s _1 and 
non-fossils with v c > 20 km s" 1 , but a gradual transition of properties from one population to the 
other. Some fossils may become more massive than v c ~ 20 km s" 1 after reionization, accrete gas 
from the IGM, and form a younger stellar population. If the dark halo circular velocity remains 
close to 20 km s _1 the young stellar population is likely to be small with respect to the old one. 
We call these galaxies "polluted fossils" because they have the same basic properties of "fossils" 
with a sub-dominant young stellar population (see RG05). 

Vice versa, some non- fossil galaxies with v c > 20 km s" 1 may loose a substantial fraction of 
their mass due to tidal interactions. If they survive the interaction, their pro perties, such as s urface 
brightness and half light radius, may be either modified or stay the same. iKravtsov et al.l (J2004) 
estimate that 10% of Milky Way dark matter satellites were at least ten times more massive at 
their formation than they are today and more were a few times more massive than they are today. 
Although their simulation does not include stars, they favor the idea that the stella r properties of 



these halos would be unchanged. Conversely, a recent work of lPeharrubia et al.l (|2008l ) looks at tidal 
stripping of dark matter and stars, achieving some success in reproducing the observed properties 
of ultra-faint dwarfs assuming that they are tidally stripped dlrrs. Using our simulations, we 
cannot make predictions of the internal properties of non-fossil dwarfs, which are too massive to be 
present in significant numbers in the small volume of our simulation. However, using our simulation 
data (in RG05), GK06 finds that about one third of Milky Way satellites may be fossils based on 
comparisons between observed and simulated Galactocentric distribution of the satellites. 
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3.1. Description of the simulation 



Given a cosmological model, simulating the formation of the first stars is a relatively well 
defined initial condition problem. However, these simple initial conditions become unrealistic as 
soon as a few stars form within a volume of several thousands of co- moving Mpc 3 . The physics 
becomes complex as competing feedback effects determine the fate of the first galaxies: radiative 
feedback regulates the formation and destruction of H2 in the intergalactic medium and in proto- 
galaxies. The background in the Lyman- Werner bands (between 11.3 eV and 13.6 eV) dissociates 
H2 through the two-step Solomon process. A few Population III stars per Mpc 3 can suppress or 
delay the formatio n of the first small mass galaxies, reducing drastically their number per unit 
co-moving volume (jHaiman et al.ll2000l ; ICiardi et al.ll2000l ; iMachacek et al.ll2000l ). 



Ionizing radiation may enhance the production of H2 (we refer to this as "positive feedback") 
by creating fr ee electrons and promoting the formation of H~, the main catalyst for the for- 



mation of H? ( Shapiro & Kane. 



1987; 



Haiman. Rees. &; Loeb 



1996; 



Rico tti. Gnedin. Shull 



2001 



Alvarez. Bromm. fe; Shapirdl2006l : ICiardi et al.ll2006l ). iRicotti et al.l (120011 ) studied the structure of 
H II regions in the early universe in a gas of primordial composition. It was found that shells of 
H2 can be continuously created in precursors around the Stromgren spheres produced by ionizing 
sources and, for a bursting mode of star formation, inside recombining H II regions. This is because 
the catalyst H~, and hence H2, is formed in regions where the gas ionization fraction is about 50%. 
The thickness of these regions depends on the density of the gas and the spectrum of the ioniz- 
ing radiation. This local "positive feedback" may be the dominant feedback, but it is difficult to 
incorporate into cosmological simulations because the implementation of spatially inhomogeneous, 
time-dependent radiative transfer is computationally expensive. 

The simulation analyzed in the present paper differs from other studies b e cause it includes 
"positive feedback" from ionizing radiation self-consistently ( Ricotti et al. 2002al lbl. I2OO8 ) (hereafter. 
RGS02a,b, RGS08). The simulation have a space resolution of 156 pc h _1 comoving (about 15 pc 
physical at z = 10) and a mass resolution of 4.93 x 10 3 M h _1 for dark matter and ~ 657 M h _1 
for the baryons. There are more than 10 5 stellar particles in our simulations at z ~ 10. The stellar 
masses are always smaller than the baryon mass resolution but can vary from ~ 0.6 M IT 1 to 
600 M Q IT 1 with a mean of 6 M h . Stellar particles do not represent real stars. In addition to 
primordial chemistry and 3D radiative transfer, the simulations include a recipe for star formation, 
metal production by SNe and metal cooling (see RGS02a for details). The code also includes 
mechanical feedback by SN explosions. However, we found that for a Salpeter IMF, the effect 
of SN e is not dominant when compared to feedback produced by ionizing radiation from massive 
stars (jRicotti et al.ll2008l ). The effect of SN explosion is somewhat model dependent and uncertain 
because it is treated using a sub-grid recipe. Hence, the simulation analyzed in this work includes 
metal pollution but not mechanical feedback by SNe. 



RGS02b have shown that the main negative feedback thought to suppress the formation of the 
first galaxies (H2 photo-dissociation) is not dominant, contrary to the results of previous studies 
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that did not include "positive feedback" from ionizing radiation. Feedback by ionizing radiation 
plays the key role, inducing a bursting star formation mode in the first dwarf galaxies. Galactic 
outflows, produced by UV photo-heating from massive stars, and H2 formation/photo-dissociation 
induce the bursting star formation mode that acts as the catalyst for H2 re-formation inside relic 
(recombining) H II regions and in the "precursors" of cosmological Stromgren spheres (i.e., positive 
feedback regions). As a result, star formation in the first galaxies is self-regulated on a cosmological 
scale - it is significantly reduced by feedback but it is not completely suppressed, even in small mass 
halos with v c ~ 5 — 10 km s . Due to the feedback-induced bursting mode of star formation in 
pre-reionization dwarfs, the cosmological H II regions that they produce remain confined in size 
and never reach the overlap phase that defines the epoch of reionization. 

The simulation data shown in this paper is the higher resolution run in RGS02b, evolved 
further to redshift z = 8.0 after the introduction of a bright source of ionizing radiation that 
completes reionization at z ~ 9 (see RG05 for details). The need for introducing a bright ionization 
source is dictated by the small volume of the simulation (1.5 3 Mpc 3 ); otherwise the volume would 
be reionized too late. The H I ionizing source removes all the remaining gas from halos with 
v c < 20 km s _1 and shuts down star formation. 

3.2. Statistical properties of simulated "fossils" vs observations 

Here we compare the RG05 predictions for the fossils of primordial galaxies to the observed 
properties (see Table H]) of the new Milky Way and M31 dwarfs. The symbols and lines in Figs.[T]l6] 
have the following meanings. All known Milky Way dSphs are shown by circles; Andromeda's dSphs 
satellites are shown by triangles; simulated fossils are shown by the small solid squares. The solid 
and open symbols refer to previously known and new dSphs, respectively. The transition between 
fossils and non-fossil galaxies is gradual. In order to illustrate the different statistical trends of 
"non-fossil" galaxies we show dwarf irregulars (dlrrs) with asterisks and the dwarf ellipticals (dE) 
as crosses, and we show the statistical trends for more luminous galaxies as thick dashed lines on 
the right side of each panel. 

Figure [T] shows how the surface brightness (top panel) and half light radius (bottom panel) of 
all known Milky Way and Andromeda satellites as a function of V-band luminosity compares to the 
simulated fossils. The surface brightness limit of the SDSS is shown by the thin solid lines in both 
panels of the figure. The new dwarfs agree with the predictions up to this threshold, suggesting the 
possible existence of an undetected population of dwarfs with below the SDSS sensitivity limit. 
The new M31 satellites have properties similar to their previously known Milky Way counterparts 
{e.g., Ursa Minor and Draco). Given the similar host masses and environments, is reasonable to 
assume a similar formation history for the halos of M31 and the Milky Way. This suggests the 
existence of an undiscovered population dwarfs orbiting M31 equivalent to the new SDSS dwarfs. 
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Fig. 1 . — An extension of Figu r e 7 from RG05 to inclu de new dwarfs i n the SPSS (Bcloku rov et al 



20071. bond: llrwin et al.l I2OO7I : Iwillman et "all l2005bl lal: Iwalsh et al 



Geha et al.l IgOOgh and recent surveys of Andromeda (jlbata et al 



2007 



2007; 



Zucker et al.1 



2006b 



Maiewski et al.l 120071 ; 



Martin et al.ll2006l ). Surface brightness and core radius are plotted vs. V-band luminosities. Small 
filled squares are the RG05 predictions, asterisks are known survivors, crosses are known polluted 
fossils, closed circles are the previously known dSph around the Milky Way, closed triangles are 
previously known dSph around M31, and open circles and triangles are new dSph around the Milky 
Way and M31 respectively. The solid lines show the SLOAN surface brightness limits and the dashed 
lines sho w the scaling relationships fo r more luminous Sc-Im galaxies (10 8 L Q ^ Lb ~ W u Lq) de- 



rived bv lKormendv &: Freeman! ([2004). 
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3.2.1. Extreme Mass-to-Light ratios 

The large mass outflows due to photo-heating by massive stars and the consequent suppression 
of star formation after an initial burst, make reionization fossils among the most dark matter 
dominated objects in the universe, with predicted M/L ratios as high as 10 4 and Ly ~ 10 3 — 10 Lq. 

Figure shows the velocity dispersion (bottom panel) and mass-to-light ratios, M/Ly (top 
panel), as a function of V-band luminosity of the new and old dwarfs from observations in com- 
parison to simulated fossils. The symbols are the same as in the previous figures. While mass data 
is available for all the previously known dw arfs, we found no published a va lues for 9 dlrr, 4 dE 



and 3 dSph (Antila, Phoenix and SagDIG) ([Mated Il998l : IStrigari et alj|2008l ). We observe a good 



agreement between the statistical properties of the new dwarf galaxies and the RG05 predictions 
for the fossils, although simulated dwarfs show a slightly larger mass-to-light ratios than observed 
ones at the low luminosity end, Ly < 10 4 L . Theoretical and observed dynamical masses are 
calculated the same way, from the velocity dispersions of stars (i.e., M = 2r 1 / 2 o" 2 /G'), and do not 
necessarily reflect the total mass of the dark halo at virialization. Indeed, the simulation provides 
some insight on why the observed value of the dynamical mass, M ~ (1 ± 5) x 10 7 M , remains 
relative constant as a function of Ly. Simulations show that in pre-reionization dwarfs, the ratio of 
the radius of the stellar spheroid to the virial radius of the dark halo decreases with increasing dark 
halo mass. The lowest mass dwarfs have stellar speroids comparable in size to their virial radii (see 
RGS08). As the halo mass and virial radius increases, the stellar spheroid becomes increasingly 
concentrated in the deepest part of the potential well, thus the ratio of the dynamical mass within 
the largest stellar orbits to total dark matter mass is reduced. This effect maintains the value of 
the dynamical mass within the stellar spheroid (measured by the velocity dispersion of the stars) 
remarkably constant even though the total mass of the halo increases. 

If dwarfs are undergoing tidal disruption (e.g., Ursa Major II), the velocity dispersions could 
be artificially inflated. However, the agreement with theory is rather good for all the new ultra- 
faint dwarfs. The data on the lowest luminosity dwarfs in our simulations are the least reliable, 
because they are very close to the resolution limits of the simulation (we resolve halos of about 
10 5 M Q with 100 particles). Let us assume we can trust the simulation and the observational data, 
even for the lowest luminosity dwarfs and the discrepancy between simulation and observation at 
the low luminosity end is real. The dynamical mass is M cx ri^cr 2 . There is good agreement 
between observations and simulations for a in faint dwarfs. Thus, it is likely that the the reason 
for disagreement in M/Ly is due to the value of ri/2, being smaller for the observed dwarfs than 
the simulated ones. This could be partly due an observational bias that selects preferentially dSphs 
with higher surface brightness and smaller r\/2 (see Figure [1]). Another explanation, is a dynamical 
effect not included in our simulation that reduces the stellar radius of dwarfs after virialization; for 
instance, tidal stripping or relaxation. A larger sample size of distant dwarfs, including kinematics 
on And XII and And XIV, both believed to be on their first approach to M31, would be useful to 
better characterize this discrepancy. 
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Fig. 2. — Mass-to-light rati o and velocity dispersion of a subset of the new dwarfs (jMartin et al 



2007 



Simon &; Gehal 120071 ) plotted with values for previously known dwarfs and RG05 pre- 



dictions. The dashed lines show the scaling relationships for more luminous Sc-Im galaxies 
(10 8 L© < L B < 10 n L Q ) derived bv lKormendv fc Freeman! ibooi l. 
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The velocity dispersions of the stars for the observed new dwarfs and the simulated fossils are 
just below 10 km s _1 for luminosities < 10 6 Ly. The scatter of the predicted and observed velocity 
dispersions as a function of dwarf luminosity suggests that dwarfs with the same luminosity may 
be found in a broad distribution of halos masses. These results are in agreement with Figure dj 
showing that the faintest primordial galaxies of a given luminosity may form in halos with a total 
mass at virialization between 10 6 — 10 7 M to a few times 10 8 Mq. This is because at these small 
masses the star formation efficiency is not necessarily proportional to the dark halo mass and there 
is a large scatter in /* at any given mass. This is due to the nature of feedback effects that is 
local and depends on the environment (e.g., positive feedback on the formation of H2). In RG08 is 
found that pre-reionization dwarfs that form in relative isolation have typically smaller value of /* 
than dwarfs of the same total mass that form in the vicinity of other luminous dwarfs (see Fig. 4 
in RGS08). 



3.2.2. Understanding the luminosity-metallicity relation 

The metallicity-luminosity relation of the observed and simulated dwarfs is shown in Figure El 
[Fe/H] is plotted against V-band luminosity in solar units. Symbols for the previously known 
dwarfs, the new, ultra-faint dwarfs, and simulated fossils are the same as in Figure HJ In this 
plot, we also color code simulated fossils according to their star formation efficiency, /*, defined as 
/* = M^/Mbar, where M* is the mass in stars and -/Vf& ar ~ M^ m /6 is the baryonic mass of the halo 
assuming cosmic baryons abundance. Red symbols show simulated dwarfs with /* < 0.003, blue 
0.003 < /* < 0.03 and green /* > 0.03. 

The new ultra-faint dwarfs do not follow a tigh t luminosity-metallicity relationship observed in 



more luminous galaxies (but see lGeha et al.1 (120091 )). This behavior is in good agreement with the 
predictions of our simulation. The recently found dwarf galaxy Segue 1 with luminosity of 34OL 
and metallicity of [Fe/H]~ —2.8 fills a gap in the luminosity -metallicity plot that was previously 



devoid of observed dwarfs, present instead in the simulation (jGeha et al.ll2009l ). However, Segue 1 
has a half light radius that is smaller than what our simulation predicts. Compared to their 
previously known counterparts, the new dwarfs have a slightly lower metallicity but much lower 
luminosities. 

There are several physical mechanisms that may produce the observed scatter in metallicities 
of dwarfs at a given constant luminosity. Here, we identify the two mechanism that are dominant 
in our simulation for primordial dwarf galaxies: 1) the large spread of star formation efficiencies 
producing a dwarf of a given luminosity is the dominant mechanism (to zero-th order approximation, 
in a closed box model, we have Z oc /*) and 2) the existence of dwarfs experiencing either a single 
or multiple episodes of star formation contribute to the metallicity spread as well. Metal pollution 
from nearby galaxies at formation might also play a role. For more massive dwarfs th at form after 



reionization, there may be different processes that dominate metal enrichment. See iTassis et al 



(2008) for a discussion. 
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Let's start from the first mechanism. Contrary to what it is usually the case for more massive 
galaxies, we find that in primordial dwarfs with masses ^5x 10 M the star formation efficiency 
/* does not monotonically increase with halo mass (i.e., /* has a large spread for a given halo 
mass or for a given mass in stars, M*, see RGS08 Figure 4 and Figure 7). The wide range in 
values results from the sensitivity of /* on a halo's environment and is due to local feedback effects 
that are of fundamental importance in determining star formation in shallow potential wells typical 
of pre-reionization dwarfs. Figure [4] (taken from Figure 7 in RGS08) shows that below a few 
10 Mq, halos with exactly the same dark mass can be either dark or luminous (depending on the 
environment). Thus, feedback from non-ionizing and ionizing UV radiation, mechanical feedback 
and chemical enrichment can produce two halos with the same dark mass and very different star 
formation efficiencies. This appears to be the main effect responsible for the observed spread in 
metallicity for a given luminosity, in the new dwarfs. In Figure El we plot the metallicity as a 
function of the mean star formation efficiency for the halo, As expected, simulated dwarfs with 
higher values of the star metallicity are the ones with the larger value of /* . 

However, this effect alone cannot account for all the observed scatter of the metallicity as 
illustrated by the color coding of simulated dwarfs in Figure El It appears that the metallicity is 
not simply proportional to /* (otherwise the boundaries between symbols of different color would 
be horizontal). Instead, for a given value of the metallicity is larger for fainter dwarfs. It 
is not too surprising that Z is not simply proportional to Even when using a very simple 
chemical evolution model, neglecting gas inflows and assuming instantaneous metal recycling, the 
mean metallicity of the stars is proportional to M*/M gas rather than /* = M*/Mb ar , where M gas 
is the initial value of the gas mass available for star formation. Thus, Z cx f*(Mb ar /M gas ). If 
feedback effects reduce the value of M gas /Mb a r below unity in the smallest and lowest luminosity 
primordial dwarfs, the metallicity of the stars will be larger for a fixed value of _/*, as observed in 
Figure [5] and Figure El The reduction of M gas /Mb ar below unity can be produced by three effects: 
the increase of the Jeans mass of the IGM over the virial mass of the halo due to reheating (see 
Figure 6 in RGS08), heating of the gas via ionizing radiation from stars within the halo, and by 
multiple episodes of star formation with a first burst that lowers M gas substantially, but does not 
produce sufficiently large values of _/* and Z when compared to subsequent bursts. 

Figure [6] shows [Fe/H] versus the surface brightness in the V-band, £y The symbols are the 
same as in Fig. [1] and the solid line shows the SDSS sensitivity limit. No trend is observed between 
metallicity and £y for the simulated dwarfs. Observed dwarfs show less scatter for Ey-metallicity 
relation than for the luminosity- metallicity relation in Figure El There is one dwarf with metallicity 
below [Fe/H ] = —2 .5: S egue- 1 that has [F e/H]= —2.8. Since the spectral synthesis method used in 



Kir by et al.l (|2008l ) and lGeha et al.1 ((2009) may not be subject to the overestimation of metallicities 
seen with measurements using the CA triplet, the lack of dwarfs with [Fe/H]< —3.0 could be a sign 
of a change in the IMF at very low [Fe/H]. 

Finally, we have examined whether there is a dependence of the metallicities on the distance of 
the galaxy from the Milky Way or Andromeda. For the new Milky Way dwarfs there is slight trend 
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of higher metallicities at smaller Galactocentric distances; however, the upward trend is dominated 
by Ursa Major II and Coma Ber., both of of which show evidence for tidal disruption. For the M31 
dwarfs, no trend was observed. 

Figures [7] show the scatter of the metallicity of the stars, o\Fe/Hh plotted against V-band 
luminosity and [Fe/H] respectively. Once more, the various point types and colors are the same 
used in Figure [3l The observational data for the new dwarfs matches the predictions, though 
Figure [7J shows a lack of low Ly objects with cj[F e /H] < 0.4. However, given the small number 
of data points available, it is not possible to rule out selections effects of small number statistics 
as an explanation. Dwarfs with low values of crpwHi tend to have higher luminosities and are 
equally likely to be a dE or dSphs. However, dwarfs with the highest cr\F e /H] are faint dSphs 
with Ly < 10 6 Lq, and occupy the lowest mass dark matter halos. It is not clear at this point 
how reliable the simulation data for cr\F e /H] is f° r dwarfs with luminosities Ly < 10 3 — 10 4 L . 
However, all the dwarfs we analyze have at least 10 stellar particles and 100 dark matter particles. 
The masses of stellar particles vary, depending on the star formation efficiency and the duration of 
the star burst. 

As with the metallicities, we looked at how the metallicity spread depends on the distance 
from the host. For both previously known and ultra-faint dSphs there is no dependence on distance 
within the virial radius of the Milky Way. There is a lack of dwarfs with high (T\F e /m beyond 
400 kpc; however, given the small number of data points and the luminosity and surface brightness 
limits of current surveys, the trend is not statistically significant. 

3.3. The Missing Galactic Satellite Problem Revisited 

In order to simulate a representative sample of the universe, the size of cosmological simulations 
must be significantly larger than the largest scale that becomes non-linear at the redshift of interest. 
At z = 0, this scale is at least 50 to 100 Mpc. Current computational resources are not able to evolve 
a cosmological simulation of this size that includes all the relevant gas and radiation physics, to 
z = 0. However, as argued in § 3, the properties of those fossil galaxies that survive tidal destruction 
change only through passive aging of their stars formed before reionization (but see Ricotti 2008), 
allowing their properties at z = to be simply related to their properties at reionization. GK06 uses 
this approximation in conjunction with high-resolution N-body simulations of the Local Group, to 
evolve a population of dwarf galaxies around a Milky Way mass halo from z = 70 to z = 0. For 
details of the simulations, see § 2 in GK06. 

GK06 define a fossil as a simulated halo which survives to z = and remains below the critical 
circular velocity of 20 km s -1 with no appreciable tidal stripping. They calculate the probability, 
Ps(v c ,max, r ), of a luminous halo with a given maximum circular velocity v c ^ max to survive from 
z = 8 (the final redshift of the RG05 simulation) to z = 0. For a given fc,max; the number of dwarfs 
at z = is N(v c ^ max , z = 8)Ps(v maX i r). The surviving halos are assigned a luminosity based on the 
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Fig. 3. — Metallicity vs. luminosity for the new and old dwarfs plotted against RG05 predictions. 
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Fig. 4. — Fraction of baryons converted into stars as function of halo mass of the galaxy at z = 10. 
Circles, from smaller to the larger, refer to galaxies with gas fractions f g = M gas /Mb < 0.1% (blue), 
0.1% < f g < 1% (cyan), 1% < f g < 10% (red) and f g > 10% (green), respectively. 
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Fig. 5. — Metallicity versus star formation efficiency /* = M^/M^ for the simulated fossils. The 
large squares show galaxies with Ly > 10 3 L , and the small squares galaxies with Ly < 10 3 L . 
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Fig. 6. — Metallicity vs. surface brightness for the new and old dwarfs. The surface brightness 
limit of the SDSS is shown by a solid vertical line. Note the predicted dwarfs with Z<-2.5 and 
surface brightnesses above Sloan detection limits. 
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Fig. 7. — Metallicity spread vs. V-band luminosity for the new and old dwarfs. Simulation data 
for the metallicity spread may be unrelayable for dwarfs with luminosities Ly < 10 3 — 10 4 Lq, due 
to the small number of stellar particles in the galaxies. 
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Fig. 8. — Luminosity function of pre-reionization fossil dwarfs predicted in GK06 (red) plotted 
with the luminosity function for the new and old Local Group dSphs. The black lines are the 
observations corrected for completeness as discussed in Section 2. Corrections to the theory account 
for an under-abundance of small halos near the hosts due to numerical effects (GK06). 
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Ly versus v c ^ max relationship from RG05. At z = 0, GK06 has a population of dwarf galaxies with 
a resolution limit of 

Vc,max — 13 km s . Unfortunately, this limit corresponds to a lower luminosity 
limit of Ly ~ 10 5 Lq, which includes Leo T and Canes Venatici I, but excludes all the other new 
ultra-faint Milky Way satellites. We do not show the prediction of the GK06 model below this 
lower luminosity limit. A new N-body method with a resolution limit of Ly ~ 10 3 L that follows 
the merger histories for the halos will be presented in Bovill & Ricotti (2008), in preparation. 

In Figure [HJ we show the cumulative luminosity function from GK06 for the Milky Way and 
M31 satellites with the addition of the new ultra- faint dSphs. The lower panel shows satellites with 
distance from their host d < 100 kpc, the middle panel d < 300 kpc and the upper panel d < 1 Mpc. 
The gray lines show the GK06 predictions, and the shaded region encompasses the error bars. Since, 
the resolution limits of GK06 causes halos with v C:max < 17 km s -1 to the preferentially destroyed 
by tidal effects, the predicted luminosity function is corrected. Both the uncorrected (lower) and 
corrected (upper) luminosity functions are plotted in the lower panel. On all three panels, the black 
histogram and the points with error bars represent the observed luminosity function of fossil dwarfs 
around the Milky Way and M31. Their numbers have been corrected for completeness as discussed 
in § 2. For the purposes of this plot, we are considering all the new dwarfs to be fossil candidates. 

For d < 100 kpc (bottom panel), there is an overabundance of observed satellites with respect 
to the simulated luminosity function from 10 5 to 10 6 L . This discrepancy is likely due to excessive 
destruction rate of satellites caused by the insufficient resolution the GK06 N-body simulations. At 
distances d < 300 kpc (middle panel), there is excellent agreement between theory and observation. 
Canes Venatici I and the new Andromeda satellites are included in the latter panel. The upper 
panel shows the luminosity function for all dwarfs within 1 Mpc of the host, including Leo T. Note, 
that GK06 assumes an isolated Milky Way type galaxy (with total mass comparable to the Local 
Group mass), while observations with d < 1 Mpc of the Milky Way include the satellite system 
around M31. For d < 1 Mpc, there is an under-abundance of observed satellites between 10 5 to 
10 6 Lq with respect to the simulation predictions. However, this is consistent with the theory since 
beyond 250 kpc dwarfs with Ly ~ 1O 5 L drop below SDSS detection limits (Koposov et al, 2007). 
Hence, the under-abundance of observed dwarfs at large distances may be due to the completeness 
limit of the survey. 



4. Discussion and Conclusions 

There are two main ideas for the origin of dSphs in the Local Group. Most importantly, these 
two ideas have very different implications for models of galaxy formation, and the minimum mass a 
dark halo needs to host a luminous galaxy. The "tidal scenario" , holds the dwarfs we see today were 
once far more massi ye, having been stripp ed of most of their dark matter during interactions with 
larger galaxies (e.g., iKravtsov et al. 1 12004 1. In this model, we would expect the halos with original 



dark matter masses below 10 8 M to be mostly dark at formation and at the modern epoch. 
The "primordial scenario", has dwarf galaxies starting with close to their current stellar mass of 
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about 10 3 — 10 6 M and, with several dark halos with mass at formation below the threshold of 
about 2 x 10 8 M hosting a luminous galaxy. Star formation in halos this small is possible only 
before reionization and is widespre ad if "positive" feedback plays a significant role in regulating 
star formation in the first galaxies (jRicotti et al.ll200ll . l2002al lbl). 



In this paper, we argue that the recent discovery of the ultra-faint dwarfs in the Milky Way 
and M31 supports the "primordial scenario". The existence of the ultra-faint dwarfs was predicted 
by simulations of the formation of the first galaxies (see RG05) and, as shown in the present work, 
the observed properties of this new population are consistent with them being the "fossils" of the 
first galaxies. 

While tidal stripping can reproduce properties of an individual galaxy, it is unable to completely 
reproduce all the trends in the ultra-faint population. This is primarily seen in the kinematics 
of the ultra-faint dwarfs. Tidal stripping predicts a steeper than observed drop in a with Ly 



(jPenarrubia et al.ll2008l ). while our simulations show primordial dwarfs which match the observed 
trends in a extremely well. It has not been shown yet that star formation in dwarf galaxies more 
massive than 10 8 — 10 9 M Q can reproduce the observed properties of ultra-faint dwarfs without 
requiring tidal stripping of stars. 

The tidal model predicts that gas rich dlrr loose their gas and transform into dSphs due to tidal 
or ram pressure interaction with a host halo. And XII, which shows a proper motion cl ose to current 



published escape velocity o f M31, may be on its first approach to the Local Group ([Martin et al 
2006; IChapman et al.l 120071 ). A similar situation exists for And XIV. With a dynamical mass of 



M ~ 3 x 1O 7 M , An d XIV has v > vM^ \ suggesting it is also just entering the Local Group 
dChapman et al.ll2007h . In the tidal model (|Maver et alibooibood ). And XII and And XIV would 
be expected to still harbor significant res ervoirs of gas, however, observations show And XIV ha s 
M HI < 3 x 1O 3 M (IChapman et al.ll2007l ) and And XII has no detected H I ((Martin et allbood ). 
If neither of these dwarfs have undergone significant tidal interactions with their hosts, as their 
velocities suggest, how did they lose their gas? Though its velocity is unknown, the recently 
discovered And XVIII (IMcConnachie et al.ll2008l ). shows the same characteristics. At a distance of 
600 kpc from M31 and 1.35 Mpc from the Milky Way, it is unlikely that And XVIII has undergone 
significant interaction with either Local Group spirals. And X VIII is classified as a dSp h with no 
detected H I and is similar to the Cetus and Tucana dwarfs (jMcConnachie et al.l 120081 ) , both of 
which are good candidate fossil galaxies (RG05). 

On the opposite end of the spectrum is the strange case of Leo T, the properties of which are 
discussed in Section 2.1.1. While, Leo T has an H I mass fraction typical o f dlrr, its other prope rties 
are indistinguishable from the other newly discovered ultra-faint dwarfs (jSimon Gehal 120071 ). all 
of which are dSph. Leo T's large dista nce from its host, H I reservoir and low probability of recent 
tidal interactions (jde Jong et al.ll2008l ) make it a good candidate for a precursor to a dSph in the 
tidal scenario. Particularly given th at Leo T's dynamical mass within the stellar spheroid is small: 
8.2 x 10 6 M ( Simon &: Gehal 120071 ) . its gas is unlikely to survive a single tidal encounter intact. 
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Therefore, Leo T may have formed at or near its current mass, and the striking similarity of Leo 
T to its ultra-faint counterparts suggests that they too could have formed as primordial dwarfs at 
their current masses. 

By our definition, pre-reionization fossils are dwarfs that form before reionization in dark halos 
with v c < v^' ~ 20 km s _1 , while non- fossils dwarfs form in halos with v c > v^ r before and after 
reionization. The value v"' ~ 20 km s~ 1 that we use to define a fossil is primarily motivated by 
fundamental differences in cooling and feedback processes that regulate star formation in these 
halos in the early Universe. This value of the circular velocity is als o very close to estimates based 



on th e suppression of star for mation in dw arfs after reionization (jGnedinl |2000| ; lOkamoto et al 



2008). However, as argued in (jRicottil 120091 ) . fossils dwarfs can have a late phase of gas accretion 
and star formation well after reionization, at redshift z < 1 — 2. Thus, a complete suppression of 
star formation after reionization is not necessarily what defines a "fossil galaxy" . 

The number of Milky Way dark satellites that have or had in the past v c > v^ r can be estimated 



using the results of published N-body simulations (see § 12. ip . We find that using the via Lactea 
N-body simu lation there are appro ximately Ndark ~ 73 ± 16 halos with v a > 20 km s -1 within the 
virial radius (jDiemand et al.ll2007l ). The new Aquarius simulations (jSpringel et al.ll2008l ). however, 
show a factor of 2.5 increase in the number of halos with v c > 20 km s _1 , i.e., Ndark ~ 182 ± 40 
dark halos. Within a distance of 200 kpc we estimate Ndark ~ 36 ± 8 for the Via Lactea and 
Ndark ~ 91 ± 20 for the Aquarius simulation. 

If the number of observed dwarf satellites within the Milky Way (after applying completeness 
corrections) is larger than Ndark we must conclude that some satellites are fossils. Twelve new 
ultra-faint dwarfs have been discovered around the Milky Way by analyzing SDSS data in a region 
that covers about 1 /5 of the sky. Applying a simple correction for the sky coverage we estimate that 
there should be about at least 85 ± 14 Milky Way satellites. However, the data becomes incomplete 
for ultra-faint dwarf that are further than about 200 kpc from the Galactic center. Comparing 
this number of luminous satellites to Ndark within 200 kpc we cannot conclusively conclude that 
some ultra-faint dwarfs are fossils because Ndark for the Aquarius simulation is comparable to the 
estimate number of luminous satellites. 



Once both sensitivity and survey area corrections are applied, iTollerud et al.l (|2008l ) estimates 
the existence of 300 to 600 luminous satellites within the viri al radius (R„j r ~ 40 kpc) of the 



Milky Way and 120 within 200 kpc. Comparing Ndark to the I Toller ud et al.l (J2008J) estimates of 
the number of luminous Milky Way satellites implies that a significant fraction of them are fossils 
(regardless if we use the Via Lactea or the Aquarius simulations estimates for Ndark ) • In Table 12.11 
we have summarized the aforementioned results. 

Another argument for the existence of fossils is provided by detailed comparison of the Galac- 
tocentric distribution of fossils in the Milky Way (GK06). Based on these comparison GK06 find 
that about 1/3 of Milky Way dwarfs may be fossils. In this paper, we show the GK06 theoretical 
results in comparison to updated observational data, including the new ultra-faint dwarfs found 



-27- 



using SDSS data, and applying completeness correction due to the limited area surveyed by the 
SDSS (about 1/5 of the sky). Assuming that the Local Group has a mass of 3 x 10 12 M , as in the 
GK06 simulation, we find that there are no "missing galactic satellites" with Ly > 10 5 L within 
the virial radius of the Milky Way. When the new dwarfs are included, the observed and predicted 
numbers of satellites agree near the Milky Way, however, for distances greater than 200 kpc, it is 
clear that there is still a 'missing' population of dwarfs. However, given that for d > 200 kpc, dwarfs 
with Ly ~ 10 5 L drop below SDSS detection limits (Koposov et al, 2007), the under-abundance 
of observed dwarfs at large distances is not surprising and likely due to the SDSS sensitivity limit. 

A final comment regards the cosmological model. The RG05 and GK06 simulations use cosmo- 
logical parameters from WMAP 1. N-body simulations show that the number N(M) of Milky Way 
dark matter satellites as a function of their mass is not overly sensitive to the cosmology, a lthough 



there are some differences on the number of the most massive satellites (Madau et al.ll2008f ). How- 



ever, N(v max ) should be sensitive to the cosmology (IZentner Sz Bullockl 12003 ). and changes of as 
and n s may affect the occupation number and Galactocentric distribution of luminous halos. The 
collapse time of small mass halos in high density regions probably dominates the 20% variations 
in o"8 between WMAP 1 and WMAP 3, limiting effects due to the cosmology near large halos. A 
decrease in luminous dwarf numbers, due to the lower erg, could be evident in the distribution of 
the lowest mass luminous halos in the voids. 

In conclusion, the number of Milky Way and M31 satellites provides an indirect test of galaxy 
formation and the importance of positive feedback in the early universe. Although the agreement of 
the SDSS and new M31 dwarfs' properties with predictions from the RG05 and GK06 simulations 
does not prove the primordial origin of the new ultra- faint dwarfs, it supports this possibility with 
quantitative data and more successfully than any other proposed model has been able to do so far. 
At the moment, we do not have an ultimate observational test that can prove a dwarf galaxy to be 
a fossil. Even a test based on measuring the SFH of the dwarf galaxies may not be discriminatory 
because, as has been recently suggested, fossil galaxies ma y have a late phase of gas accretion and 



star formation at z < 1 — 2, during the last 9 — 10 Gyrs (lRicottill2009l ). The distinction between 
fossils and non-fossils galaxies thus is quite tenuous and linked to our poor understanding of star 
formation and feedback in dwarf galaxies. Arguments based on counting the number of dwarfs in 
the Local Universe probably provide the most solid argument to prove or disprove the existence of 
fossil galaxies. In the future, a possible test may be provided by deep surveys looking for ultra-faint 
or dark galaxies in the local voids. Some fossil dwarfs should be present in the voids if they formed 
in large numbers before reionization. 
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Table 2. Summary of New Dwarfs Properties 
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Table 3. Summary of New Dwarfs Properties cont. 
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